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Summary. - A formal supersymmetric quantum mechanics (SUSY QM) 
procedure for any linear homogeneous second-order differential equation is 
briefly sketched up and applied to a simple exactly solvable case. 



PACS 11.30.Pb - Supersymmetry. 



It is quite well known how to get Riccati equations from the general form 
of homogeneous linear differential equation of the second order ]l| 

A(x)^ + B(x)^ + C(x)u = . (1) 
ax z ax 

It requires the transformation 

1 Ida 

where R{x) is an arbitrary function. Using u = Ruy and u = Ru{y + 
Ry 2 + ^y) in Eq. (£]) leads to 

dy ( R' B\ „ 2 C 

Moreover, since R(x) is an arbitrary function, one can choose it such 
that the coefficient of y be zero, i.e., one can end up with the simplified 
Riccati equation 
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dx °^ ARq ' ^ ^ 

where Ro{x) = exp[— J X (B/A)dx'}. 

Our point here is to notice that Eq. (Q) allows a direct connection with 
various SUSY QM schemes || by means of the change of independent vari- 
able z(x) = J x Rodx . This leads to a Riccati equation of SUSY QM type 
at zero factorization energy along the curve z(x) 

£+,"-*(,), (5, 

where V\{z) = — A ^^}^ can be interpreted as a Schrodinger 'potential'. 
The SUSY partner Riccati equation will be 

-g+.f-WM. (6) 

where the partner 'potential' V^z) is Darboux 'isospectral' with respect to 
Vi, i.e. 

V 2 = V 1 -2D 2 [ln^(z))] , (7) 

where D = 4- and ip is a particular solution of D 2 tp — V\ij) = 0. 

In addition, one can think of the SUSY QM scheme based on the general 
Riccati solution, as first tackled by Mielnik for the harmonic oscillator case 
||. In the latter approach, one gets a one-parameter family of Darboux 
strictly isospectral 'potentials' given by 

Ui (z; X) = Vi - 2D 2 [\n(I(z) + A)] , (8) 

where I(z) = J z ip 2 (t)dt. Moreover, 

is the modulated Schrodinger zero mode implied by this scheme in which 
the Riccati integration constant A is kept as a free parameter. In general, 
to get continuous solutions ip(z,\), one should take care of not having a 
zero denominator. This leads to conditions on the possible values of A. 
On the other hand, if one works with polynomial solutions there will be 
singularities in the log derivative of i/j. SUSY partner potentials based on 
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the n'th excited state of a Schrodinger discrete spectrum problem split in 
n + 1 branches separated by the n singularities of the log derivative [Q] . 

Even more general schemes, such as higher-order intertwinings can be 
applied leading to a rich class of 'isospectral' solutions. 

The main drawback of the formal scheme as presented in Eqs. (0)^(0) is 
that it is not at all easy to implement in practice because onto the z axis 
one does not get an easily solvable eigenvalue problem. Let us take as an 
example the Hermite equation 

u - 2xu + 2nu = . (10) 

For this case, A(x) = 1, B{x) = —2x, and C{x) = 2n. Thus, Rq{x) = e x 
and z(x) = J x e* dt. The potential V\(z) = —2n/e 2z does not lead to an 
exactly solvable problem. 

We provide in the following an exactly solvable case focusing on the 
scheme based on the general Riccati solution. Let us take the initial equation 
of the form 

D 2 x u+ (l/x)D x u + u = (11) 

Then V± = —z 2 that corresponds to an eigenvalue problem in terms of the 
Bessel functions, namely ip(z) = C\yfzJ\ (z 2 /2) + Ci^fzYx (z 2 /2). In Fig. 1 
we present a common three-dimensional plot of the formulas (|8|) and (|9|) for 
the pair of superposition constants C\ = and Ci = 1. On the other hand, 
in Fig. 2 we present a full axis one-dimensional plot for A = 0.2 and C\ = 1 
and C 2 = 0. 

One can also ask what is the equation in the initial x-axis corresponding 
to the potential Vi(z, A). One can easily show that such an equation has the 
following form 

Dl^(z(x), A) + ~D x ,/>(z(x), A) - ^^^(z(x), A) = . (12) 
For this case, the connection between the two axes is given by z = lnx. 
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Figure 1: The potential V\{z, A) and the eigenfunctions i/)(z,\) correspond- 
ing to ip(z) ~ yfz Yi (z 2 /2) for A G (0,30). The dependence on A is almost 
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flat in the range of large A. 
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X=0.2 




Figure 2: The eigenfunction tjj(z,X) corresponding to i/j(z) ~ • v /zJi(z 2 /2) 
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and the potential V\(z, A) for A = 0.2. The wave function seems to approach 
a finite value for z = 0, which can be explained by the presence of the 
opposite divergences in the potential V\(z, A) at z = 0. 
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